Evolution of field spiral galaxies up to redshifts z = 1 
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ABSTRACT 

We have gained intermediate-resolution spectroscopy with the FORS instru- 
ments of the Very Large Telescope and high-resolution imaging with the Ad- 
vanced Camera for Surveys aboard HST of a sample of 220 distant field spiral 
galaxies within the FORS Deep Field and William Herschel Deep Field. Spa- 
tially resolved rotation curves were extracted and fitted with synthetic velocity 
fields that take into account all geometric and observational effects, like blurring 
due to the slit width and seeing influence. Using these fits, the maximum rota- 
tion velocity V max could be determined for 124 galaxies that cover the redshift 
range 0.1 < z < 1.0 and comprise a variety of morphologies from early-type 
spirals to very late-types and irregulars. The luminosity-rotation velocity dis- 
tribution of this sample, which represents an average look-back time of ~5 Gyr, 
is offset from the Tully-Fisher relation (TFR) of local low-mass spirals, whereas 
the distant high-mass spirals are compatible with the local TFR. Taking the 
magnitude-limited character of our sample into account, we show that the slope 
of the local and the intermediate- z TFR would be in compliance if its scatter 
decreased by more than a factor of 3 between z ~ 0.5 and z ~ 0. Accepting 



1 Based on observations collected at the European Southern Observatory, Cerro Paranal, Chile (ESO 
Nos. 65.O-0049, 66.A-0547, 68.A-0013, 69.B-0278B and 70.B-0251A) and observations with the NASA/ESA 
Hubble Space Telescope, PID 9502 and 9908. 
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this large evolution of the TFR scatter, we hence find no strong evidence for a 
mass- or luminosity-dependent evolution of disk galaxies. On the other hand, 
we derive stellar M/L ratios that indicate a luminosity-dependent evolution in 
the sense that distant low-luminosity disks have much lower M/L ratios than 
their local counterparts, while high-luminosity disks barely evolved in Mj L over 
the covered redshift range. This could be the manifestation of the "downsizing" 
effect, i.e. the succesive shift of the peak of star formation from high-mass to 
low-mass galaxies towards lower redshifts. This trend might be canceled out in 
the TF diagram due to the simultaneous evolution of multiple parameters. We 
also estimate the ratios between stellar and total masses, finding that these re- 
mained constant since z — 1, as would be expected in the context of hierarchically 
growing structure. 

Subject headings: galaxies: evolution — galaxies: kinematics and dynamics - 
galaxies: spiral 



Introduction 



The concept of Cold Dark Matter (CDM) and its prediction of hierarchical structure 
growth has become an astrophysic al paradigm in the l ast decade. Observations of the Cos- 



mic Microwave Ba ckground (e.g. ISpergel et. all 120031 ) or the Large Scale Structure (e.g. 



Bahcall et al.l Il999l ) strongly support the idea that the vast majority of matter is non- 
baryonic, non-luminous and interacts only gravitationally. In this picture, gas and stars are 
embedded in Dark Matter halos, and low-mass systems were the first virialized structures in 
the early cosmic stages, followed by an epoch of accretion and merger events during which 
larger systems were successively built up. Disks are destroyed during the frequent mergers 
of similar-mass ga laxies at z > 1 but can regrow via accretion events at lower redshifts (e.g. 
Abadi et alil2003h . 



An important tool for the quantitative test of the predictions from numerical simu- 
lations based on the hierarchical scenario are scaling relations that link galaxy parame- 
ters. Within the last years, several studies focused on the evolution of such relations up 
to redshifts z ~ 1 - corresponding to more than half the age of the universe - or even 
beyond. Many analyses utilized the Tully-Fisher Relation ( TFR) between the maximum ro- 
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Studies focusing on the stellar mass TFR consisten t ly have found no e vidence for evo- 



lution over the past ~7 Gyr (e.g. IConselice et al.ll2005l ; iKassin et al.l 120071 ). On the other 
hand, analyses exploiting the rest-frame 5-band TFR — which is less sensitive to the stel- 
lar mass than to the fraction of young stars — yielded quite discrepant results. This is 
somewhat suprising since most recent TFR studies comprise ~ 100 objects or more. Some 



authors derived much higher luminosities of spiral galaxies in the past (e.g. iBohm et al. 



in luminosity (e.g. IVogt 



2001 



Flores et al. 



2004 ; iBamford et aL 120 06; Chi u et al. 2007) whil e others find only a very modest evolution 



20061 ) . There also still is debate on whether the 



i?-band TFR changes its slope with redshift whic h woul d imply an evolution depending on 
Vmax and hence total galaxy mass. In lBohm et al.l (120041 . hereafter B04), we found that low- 
mass spirals at (z) w 0.5 were brighter by up to several magnitudes than in the l ocal universe, 
while high-mass spirals barely evolved in Mb at given V m3iX . To the contrary, IWeiner et al. 



(120061 ) observe a stronger brightening in the high-mass than in the low-mass regime to- 
wards larger redshifts. From the theoretical side, the situation is similarly unclear: in some 
simulations, the 5-band TF R slope remains constant with lookback time and only the magn i- 
tude zero point changes (e.g. ISteinmetz fc Navarro! 1 19991 ; iPortinari fc Sommer-Larsenl 12 007) 



while also an evolutio n in slope is predicted by other authors (e.g. iBoissier fc Prantzos 
Ferreras fc SilklhoOlh . 



2001 



A growing number of observational studies based on other methods than the TFR 
have pointed to wards a mass-depend ent evolution of distant ga laxies, either in terms o f 
their colors (e.g. Kodama et al. 2004]). mass-to- light ratios (e.g. van der Wei et al. 2005 ). 



or average stellar ages (e.g. iFerreras et al.l 120041 ). These results indicate that the global 



stellar populations of distant high-mass galaxies are on average older than those of distant 
low-m ass galaxies, similar to what has been found in the local universe (e.g. iBell fc de Jong 
20001 ) . Since small galaxies are understood as ancient building blocks within the framework 
of hierarchical structure growth, and in turn should have an older stellar content than larger 
systems formed more recently, these observations are at variance with the straightforward 
expectation. They can be interpreted such that the peak of star formation shifts from 
high-mass to low -mass galaxi e s with growing cosmic age, a phenomenon that was termed 
"downsizing" by ICowie et al.l (119961 ). It is however a non-trivial question what kind of 
imprint "downsizing" has on the TFR since the luminosity- V m3iX diagram most likely probes 
the evolution of various galaxy properties simultaneously, including stellar M/L ratio, gas 
mass fraction, dust content etc. 

In this paper, we report on an extensive observational study of spiral galaxy evolution 
over the last 7 Gyr. The new sample holds 124 intermediate- z disk galaxies with reliable 
Vmax measurements, roughly doubling the size of our previous data set from B04. The 
analysis given therein was limited to ground-based imaging, while we here can rely on our 
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high-resolution imaging performed with the Advanced Camera for Surveys of the HST. This 
allows a much more accurate determination of disk inclination angles, scale lengths etc., and 
in turn a more robust estimate of the maximum rotation velocities. Extending the approch 
used in B04, we will analyse the evolution not only of the TFR but also of the stellar mass-to- 
light ratios and the stellar mass fractions. We stress that this is the first re f ereed publication 
on the total sample; only some brief excerpts were shown in IZiegler et al.l (120061 ). 

The structure of this paper is as follows. In §2, we describe the target selection, spec- 
troscopic data reduction and redshift determination. §3 briefly introduces the HST /ACS 
imaging and morphological analysis. Details on the derivation of the intrinsic maximum 
rotation velocities and luminosities are given in §4. We present and discuss our results in §5, 
followed by a short summary in §6. Throughout this article, we will assume a flat cosmology 
with Vt x = 0.7, Vt m = 0.3 and H = TOkms" 1 Mpc" 1 . 



2. VLT spectroscopy 



We selec ted our targets ut ilizing the multi-band photometric surveys in t he FORS Deep 



Field (FDF, Heidt et al. 



20031 ) and William Herschel Deep Field (WHDF, Metcalfe et al. 



200ll ). Note that the sample used in B04 comprised only FDF data. The following criteria 
were applied: a) total apparent brightness R < 23 m , b) spectrophotometr ic type "later" than 
Sa, based on a photometric redshift catalogue in the case of the FDF (see iBender et al.ll200ll ) 
or color-color diagrams for the WHDF objects (n o photometric r edshif ts were available for 
these), adopting the evolutionary tracks given in iMetcalfe et al.l (120011 ). c) disk inclination 
angle % > 40° and misalignment angle between apparent major axis and slit direction 5 < 15°. 
The two latter constraints were chosen to limit the geometric distortions of the observed 
rotation curves. Apart from these limits, our selection was morphologically "blind". No 
selection on emission line strength was performed. 

Using the FORS 1 & 2 instruments of the VLT in multi-object spectroscopy (MOS) 
mode, we took spectra of a total of 220 late-type galaxies between September 2000 and 
October 2002. In t he MOS masks, we also included a total of ~40 early-type galaxies for a 
different study (see lZiegler et al.ll2005l ). In the case that a target initially selected as an E/S0 
galaxy turned out to have the spectrum of an early-type spiral (Sa/Sab), it was included in 
the disk galaxy sample a posteriori. Such a strategy helps to avoid selection biases against 
high M/L ratios. 

We used a fixed slit width of 1.0 arcsec which resulted in spectral resolutions of R ~ 1200 
for the FDF observations, where the grism 600R was used, and R ~ 1000 in the WHDF 
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with the grism 600RI, respectively. Each setup was exposed for a total of 2.5 hours under 
seeing conditions between 0.43arcsec and 0.92arcsec FWHM, with a median of 0.76arcsec. 
The spatial sampling in the final spectra was 0.2 arcsec/pixel (FDF) and 0.25 arcsec/pixel 
(WHDF, after the FORS 2 CCD upgrade). 

The data reduction was conducted on single extracted spectra. All reduction steps were 
performed on the individual exposures before they were finally combined. The typical rms of 
the dispersion relation fitted for wavelength calibration was 0.04 A. For 202 of the galaxies, 
a redshift determination and spectral classification was feasible. The objects range from 
z = 0.03 to z — 1.49 with a median (z) = 0.43 (the redshift distribution of the final sample 
entering the Tully-Fisher analysis is given in §4). 



3. HST/ACS imaging 

The HST/ACS observations were carried out during cycles 11 and 12. To cover the 
6x6 arcmin 2 sky areas of the FORS and William Herschel Deep Fields, 2x2 mosaics were 
taken with the Wide Field camera that has a field-of-view of ~ 200 x 200 arcsec 2 and a pixel 
scale of 0.05 arcsec. 4 visits each with a total exposure time of 2360 s (FDF) and 2250 s 
(WHDF) through the F814W filter were used per quadrant. Each visit was split into two 
exposures. We kept the resulting frames from the ACS pipeline reduction (including bias 
subtraction, flatfielding and distortion correction) and used a filtering algorithm to combine 
the two exposures of each visit and remove the cosmics. 

The structur al parameters of the galaxies were derived with the GALFIT package 



( iPeng et al.ll2002l ). For convolution of the model profiles, a mean Point Spread Function 
was constructed using ~ 20 unsaturated stars with I^u < 23 which were normalized to the 
same central flux and median-averaged. An exponential profile was used to fit the galaxies' 
disk components, while a Sersic profile was taken to model an additional bulge, where de- 
tectable. Note that for the analysis presented here, the most important parameters are the 
inclination, position angle and scale length of the disk. The best-fitting bulge parameters 
were only used for a raw morphological classification. 

We found significant bulge components in only 46 out of 124 galaxies that were reliable 
for V max determination (see next section), spanning the range 0.01 < B/T < 0.53 (median 
(B/T) = 0.15). In terms of a visual classification within the Hubble scheme, our sample 
comprises all types from Sa spirals to irregulars. 
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4. Tully Fisher Parameters 

To extract rotation curves (RCs) from the 2D spectra, Gaussian fits were applied to the 
emission lines stepwise perpendicular to the direction of dispersion. An averaging boxcar of 
3 pixels widths, corresponding to 0.6arcsec, was used to increase the S/N. For very weak 
lines, this boxcar size was increased to 5 pixels. In cases of multiple usable emission lines 
for a given object, the RC with the largest spatial extent and highest degree of symmetry 
was used as reference. RCs from different emission lines agreed within the errors (at least in 
terms of the derived maximum rotation velocities) for the majority of the objects. Most of 
the "reference" RCs were extracted from the [0 II] or [0 III] lines, a smaller fraction from 
H a or Hg (for a vast number of example rotation curves, please see B04). 

Due to the small apparent sizes of the spirals in our sample, the observed rotation curves 
were heavilly blurred. At redshifts of z ~ 0.5, the apparent scale lengths are of the same 
order as the slit width (one arcsec for the entire data set) and the seeing disk. The maximum 
rotation velocity hence cannot be derived "straightforward" from the outer regions of the 
observed rotation curves. 

To account for the blurring effects, we generated synthetic rotation velocity fields ac- 
cording to the measured disk inclination, position angle and scale length of a given object. 
We assumed an intrinsic linear rise of the rotation velocity at small radii which turns over into 
a regime of constant rotation velocity at large radii, with the turnover radius computed from 
the scale length of the emitting gas (for details on this prescription as well as tests of various 
rotation curve shapes, we again refer to B04). Using the position angle and disk inclination 
derived with GALFIT, a 2D rotation velocity field then was generated. After weighting 
with the luminosity profile and convolution with the Point Spread Function, a "stripe" was 
extracted from the velocity field that corresponded to the orientation and width of the given 
slitet used for spectroscopy. This simulation yielded a synthetic rotation curve which was 
fitted to the corresponding observed curve to derive the intrinsic value of V max . 

Following this approach, V max values could be determined for 124 spirals in our sample 
(hereafter, we will refer to this kinematic data set as the TF sample). 78 galaxies were 
rejected, partly due to disturbed kinematics, which could introduce systematic errors in 
the TF analysis. Other objects were not suitable due to low S/N, solid-body rotation or 
the lack of significant rotation within the measurement errors. The TF sample galaxies 
span the redshift range 0.05 < z < 0.97 (median 0.45, corresponding to a lookback time of 
4.7 Gyr) and comprise 19 objects with spectrophotometric type Sab, 65 Sc spirals and 40 
Sdm galaxies. V m&x values fall in the range 22km/s < Kiax < 450km/s (median (V ms , x ) = 
135km/s). 
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The computation of the luminosities benefitted strongly from the multi-band imaging 
of our targets. Depending on the redshift of a given object, the filter which best matched the 
rest-frame S-band was used to compute the absolute magnitude M B . This way, systematic 
/^-correction errors due to dependence on the SED are very small (<0.1 mag for all types 
and re dshifts in our data set). Intrinsic absorption was corrected following iTully fc Fouque 
(119851 ). In this approach, the amount of absorption only depends on the disk inclination. For 
testing purpose s, we altern a tively applied the inclination- and mass-dependent correction 
factors given by lTully et al.l (119981 ). However, this did not change the results presented in the 
following. It only is crucial that any distant sample is corrected for internal a bsorption in the 
same way as the local sample it is compared to. Since we use the data from lPierce &: Tully 
(119921 ) as local reference here, who adopted the Tully & Fouque approach, we will keep this 
convention in the following. The TF galaxies have absorption-corrected absolute magnitudes 
between M B = -16.8 and M B = -22.7 (median (M B ) = -20.3). The 78 galaxies with RCs 
that were not reliable for a V max determination have slightly lower luminosities, with a range 
-22.4 < M B < -14.8 and a median value of (M B ) = -19.8. 



In Fig. [T], we compare our sample to the local Tully-Fisher relation from lPierce &: Tully 



( 119921 . PT92 hereafter). On average, the FDF & WHDF spirals are more luminous than their 
local counterparts by (AM B ) = — 0.84 m , which could be due to younger stellar populations, 
i.e. lower M/L ratios of the distant galaxies. This interpretation is supported by the fact that 
the amount of the distant galaxies' brightening AM B increases with redshift. A linear x 2_ fit 
yields AM B = (—1.22 ± 0.40) z — (0.02 ± 0.20) [mag]. Here, we took the error propagation 
of Knax errors as well as the typical errors for the l ocal spirals into account. Th is luminosity 
evolution is in agreement with the findi ngs by , e.g.. iMilvang- Jensen et al.l ( 120031 ). but exceeds 
the very modest TF offsets derived by I Vogt I (120011 ). 



5. Discussion 

The luminosity-rotation velocity distribution of our sample indicates that the offsets 
from the local Tully-Fisher Relation change not only with redshift (see above) but also 
with mass: while the high-mass spirals are in relatively good agreement with the local 
TFR, the low-mass spirals are overluminous by up to several magnitudes at given V max 
(cf. Fig. []JJ. Using a 100 iteration bootstrap bisector fit, we find an intermediate-redshift 
TFR of M B = -(4.27±0.30) log V max - (11. 18± 0.65), which is a significantly shallower slope 



(a = —4.27) than locally, where a = —7.48 is observed (IPierce fc Tullylll992l ); we emphasize 
that this published value is in good agreement with a bisector fit to the local sample which 
yields a = —7.57 ± 0.38. Note that the bisector fitting method - which is a combination of 
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1.5 2 2.5 



lo g V max [km/s] 

Fig. 1.— The distant FDF & WHDF galaxies at (z) = 0.45 compared to the local Tully- 
Fisher relation as given by Pierce & Tully (1992, solid line; the dashed lines denote the 3 a 
limits). The distant spirals are systematically overluminous for their values of Vm ax . Typical 
error bars for the high-mass and low-mass regime are indicated in the upper right and lower 
left corner, respectively. 
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a "forward" and an "inverse" TF fit - can be only we a kly a ffected by potentially correlated 



errors in AM B and V max suspected by iBamford et al.l (120061 ) 



We reported earlier on a potential mass-dependent luminosity evolution with an anal- 
ysis that was limit ed to the FDF sub-sample and ground-based structural parameters 



(IZiegler et al.ll2002l . B04). Therein, we ruled out a variety of systematic errors that poten- 
tially could bias the observed distant TFR slope. Firstly, we tested whether tidally induced 
star formation in close galaxy pairs affects our sample, but found that the TF distributions 
of pair candidates and isolated galaxies are consistent. Secondly, our results were robust 
against the use of different prescriptions for the intrinsic absorption correction. Thirdly, we 
sub-divided the FDF sample according to the rotation curve quality in terms of radial extent 
and symmetry. A re-analysis using only high-quality rotation curves confirmed the shallower 
distant TF slope, showing that our findings are not induced by perturbated or truncated 
rotation curves. 

We also tested whether a slope change could be mimicked by an incompleteness effect 
arising from the apparent magnitude limit in our target selection. Towards higher redshifts, 
such a limit corresponds to higher luminosities and, in turn, higher masses. A fraction 
of the low-luminosity, low-mass (slowly rotating) spiral population is therefore missed in 
the selection process, while the low-mass galaxies that are selected might preferentially be 
located at the high-luminosity side of the relation. 

To test this effect on the new, full sample, we split it into three equally-sized redshift bins 
corresponding to median redshifts (z) = 0.23, (z) = 0.45 and (z) = 0.73, respectively (see 
Fig. [2]). In all redshifts bins, at given V max the distant galaxy sample shows an overluminosity 
of the low-mass galaxies compared to the local TFR, while the distributions are similar at 
the high-mass end. Using bootstrap bisector fits, we find significant deviations from the 
local TFR slope in all redshift bins: a = -4.18 ± 0.35 at (z) = 0.23, a = -2.77 ± 0.56 
at (z) = 0.45 , and a = —3.00 ± 0.50 at (z) = 0.73, respectively. Also shown in Fig. [2] 
are the absolute magnitudes to which our selection criterion R < 23 m corresponds at the 
median redshift of each sub-sample (dot-dashed lines). To derive these limits, w e computed 



the fc-correction using a synthetic spectrum of type Sc by iMoller et al.l (120011 ). Note that 
the fc-korrection for transformation from R Q ^ S to -B res t is only weakly depending on SED type 
at redshifts 0.3 < z < 0.8. To achieve consistency, we corrected for intrinsic absorption at 
a disk inclination angle of i = 60°, which is the average of our data set. It is evident from 
Fig. [2] that the magnitude limit of our survey affects the covered luminosity ranges in the 
intermediate- and high-redshift bin. In the low-redshift bin ((z) = 0.23), the lack of galaxies 
with Mb > — 18 seems less likely to be induced by the R < 23 m selection. In the following, 
we aim at a quantitative estimate of the impact of sample incompleteness on the TF analysis. 
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1.5 2 2.5 



lo g V max [km/s] 

Fig. 2. — The FDF & WHDF galaxies compared to the local TFR as given by Pierce & 
Tully (1992, solid line; the dashed lines denote the 3 a limits). The distant sample is divided 
into three equally-sized sets with median redshifts (z) = 0.23 (solid circles), (z) = 0.45 
(open triangles) and (z) = 0.73 (open circles), respectively. The horizontal dot-dashed lines 
depict the limit in £>-band absolute magnitude corresponding to our selection on apparent 
brightness R < 23 m . An SED corresponding to Hubble type Sc has been assumed here. 
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The key factor for the strength of this selection effect is the scatter of the TFR, which 
has a value of Q"r = 0.41 m in the .B-band of the local PT92 sample. In a previous analysis 
( IZiegler et al.l 120021 ). we assumed that this scatter increases by a factor of 1.5 between z ~ 
and z = 0.5 due to, e.g., a broader distribution in star formation rates at earlier comic 
stages. Here, we will use a different approach by testing how strongly the TF scatter would 
have to evolve over the past ~5 Gyr for the V max -dependent TF offsets to be attributed to an 
incompleteness effect. 

We hence assumed that the slope of the TFR remained constant over the redshift range 
under scrutiny here, an d that its scatter has bee n larger at earlier cosmic times. Similar to the 
technique described by iGiovanelli et al.l (119971 ). a Schechter luminosity function (LF) form 
was fitted to the observed luminosity distribution, with the characteristic luminosity M* and 
the space density 0* as free parameters. At the faint end of the luminosity function, where 
our sample is incomplete, a slope of a = —1.2 was adop ted, which is a typical value found in 



studi es of the S-band LF at intermediate redshifts (e.g. iGabasch et al.ll2004l ; iGiallongo et al. 



20051 ). We could not determine the faint end slope directly from the data since the luminosity 
distribution of our sample peaks at Mb ~ —20 and falls off towards fainter magnitudes due 
to incompleteness. The ratio between the observed LF and the best-fit Schechter LF was 
computed and expressed as an incompleteness function y = /(Me) that ranged from y = 1 
in the case of 100% completeness to y — in the case of 0% completeness. On the basis 
of the local TFR, y = /(Mb) was then converted to y = f(V mSuX ). This was done because 
the distribution in V max is much less affected by the apparent magnitude limit than the 
distribution in M B , hence the incompleteness bias was computed as a function of V max , not 
clS cL function of Mb- 

To derive the impact of the sample incompleteness, an unbiased TFR of the form 
M B = alogVmax + b was assumed with a fixed slope a = —7.48. The TFR intercept b 
was determined implicitely from the best-fit value of the charateristic luminosity M* . We 
stress again that the purpose of this approach was only to test whether a time-independent 
TF slope could be consistent with our data set. For each object in the TF sample with 
an observed maximum rotation velocity V max , the incompleteness function y = /(Vmax) was 
taken as the probability that the given galaxy would enter the observed TF sample with 
an observed, biased a bsolute magnitude Mp = alogVmax + b ± a{V mij _ x ). As observed in 
the local universe (e.g. IGiovanelli et al.lll997l ). the scatter cr(V max ) was assumed to be larger 
for lower mass galaxies (at, e.g., V max = lOOkm/s, the scatter was taken to be a factor 
of 1.6 larger than at V ma _ x = 300km/s) with an average value satisfying (a(y max )) = o B . 
This computation of M B including a TF scatter was iterated 800 times for every object, 
yielding a simulated, biased TF relation. The difference between the absolute magnitude 
Mb = /(Vmax) of the un-biased TFR and the simulated absolute magnitude M B = f(V max ) 
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of the biased TFR was taken as the correction factor to de-bias our TF sample. Note that the 
de-biased magnitudes are fainter than the observed values and that the TFR becomes 
steeper after de-biasing. 

These computations were performed in three variants using i) the local scatter as, ii) a 
2 x local scatter and iii) a 3 x local scatter for the intermediate- z TFR. The corresponding 
de-biased TF samples were fitted with bootstrap bisector fits. To account for the fact that 
the distant galaxies cover a range in redshift, corresponding to various lookback times, the 
multiplicating factor for introducing the scatter in the simulated distant TFR was not kept 
constant for the whole sample but computed individually for each object in such a way that 
the average scatter was either equal to o~b, 2o~b, or 3 as- 

Using this approach, we found that the un-biased distant TFR slope would be a = 
—4.88 ±0.29 if the distant TF scatter was the same as locally. If the scatter doubles between 
z = and z ~ 0.5, the incompleteness effect becomes stronger and the de-biased distant 
TFR slope would be a = —6.09 ±0.28. Assuming a three times larger scatter than locally we 
found a de-biased distant slope of a = —7.64 ±0.27. However, the local TFR is also affected 
by a magnitude limit. The de-biased local TF bisector fit slope we find is a = —8.02 ± 0.41 
(of course, the unchanged local TF scatter of o~b = 0.41 m was used here). We hence conclude 
that the distant and local TFR would be consistent in terms of an incompleteness effect only 
if the TF scatter evolved by more than a factor of 3 over the past ~5 Gyr. This would be 
in agreement with the observed distant scatter which is as ~ 1.2 mag. On the other hand, 
this might be an overestimate since the free-fit scatter (i.e., with a slope a = —4.27, see 
beginning of this section) is only ~ 0.9 mag. 

These results imply that the differences between the distant and local TF distributions 
can at least in part be attributed to an incompleteness effect. However, the TFR directly 
traces only the overall luminosities of a galaxy sample, not the properties of the stellar 
populations. To perform a more detailed analysis, we converted the rest -frame, absorption- 
corrected B — R colo rs of the galaxies in our sample into fT-band M/L ratios following 



Bell fc de Jongi (120011 ). In Fig. [3], we show the Mj L K ratios versus measured K-band absolute 
magnitudes of the distant TF spirals in comparison to local disk galaxies. We have split the 
distant galaxies into two redshift bins covering z < 0.45 (filled symbols) and z > 0.45 (open 
symbols), respectively. It is evident that the distant TF galaxies on the average have lower 
mass-to-light ratios at a given luminosity than their present-day counterparts (shaded area; 
note that several of the distant galaxies are at the upper limit of the local distribution, 
indicating that our selection did not miss disks with high Mj L ratios) . 

Moreover, there is a dependency on luminosity: low-luminosity galaxies seem to evolve 
stronger in Mj Lk than high-luminosity galaxies. One possible explanation for this could be 
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Fig. 3. — Stellar mass-to-light ratios of our distant TF sample galaxies at 0.1 < z < 0.45 
(filled circles) and 0.45 < z < 1.0 (open circles), compared to the parameter range covered by 
present-day spirals (shaded area) from lBell &: de Jong! (120011 ) . The data indicate a stronger 
evolution in M/L for low-luminosity galaxies. 
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Fig. 4. — Stellar mass-to-light ratios of the distant galaxies which were not usable for a 
determination of V m3X , splitted into two redshift bins for 0.1 < z < 0.45 (filled symbols) 
and 0.45 < z < 1.0 (open symbo l s), re spectively. The parameter range covered by present- 
day spirals from iBell fc de Jong! (120011 ) is indicated by the shaded area. Similar as for the 
TF objects (see Fig. the data indicate a stronger evolution in M/L for low-luminosity 
galaxies. 
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a larger fraction of young stars in the distant low-luminosity galaxies than high-luminosity 
galaxies. Since the i^-band luminosity is a good tracer of stellar mass, this would imply 
higher luminosity-weighted stellar ages towards higher stellar masses. 

For comparison, we also computed the luminosities and mass-to-light ratios of the disk 
galaxies which were not included in the TF sample due to kinematic disturbancies, low S/N 
ratios of their emission lines or lack of significant rotation. These galaxies are shown in 
Fig. HI divided into two z-bins as in Fig. [3J The non-TF galaxies show a similar trend as the 
distant TF spirals, i.e. a stronger evolution of the M/L ratios in the low-luminosity regime 
with respect to the high-luminosity regime. Note that there is no evidence for very blue, 
luminous disk galaxies that have been excluded from the TF analysis (in principle, such a 
population of luminous disks with blue colors and irregular kinematics could be expected 
within the framework of the hierarchical scenario, e.g. due to high-mass spirals which have 
recently undergone merger events). 

Could the observed evolution of the M/L ratios indicate a difference in stellar ages 
between high- and low-luminosity disks in the distant universe? Indeed, this interpretation 
gains support from an analysis of the broad-band colors of galaxies at z > 0.5 from our data 
set with single-zone models of chemical enrichment, which yielded evidence for a dependence 
of the mean stellar ages on V max (and hence to tal mass) : the high- m ass galaxies have older 
stellar populations than the low- mass ones (see iFerreras et al.l 12004 ). This indication for an 
anti-hierarchical evolution of the baryonic compon ent of galaxies ( "downsizing" ) has also 



been found in other studies of distant galaxies (e.g. iKodama et al.l 12004 ; Ivan der Wei et al. 
20051 ). It is possible that the various evolutionary effects coming into play between the earlier 
and the local universe - the evolution of the stellar populations, gas mass fractions, dust 
content etc. - balance each other in such a way that the downsizing phenomenon is not 
reflected in a significant differential evolution of the Tully-Fisher Relation. 

The lack of hierarchical evolution in terms of the stellar populations raises the question 
whether a hierachical evolution of the Dark Matter Halos — the most fundamental prediction 
of (A)CDM cosmo logy — can be establ ished with the data. To test this, we will use a 



similar approach as lConselice et al.l (120051 ) by focussing on the evolution of the ratio between 
stellar and total mass since redshift z — 1. Based on the K-bsaid mass-to-light ratios, 
we transformed the absolute K magnitudes into stellar masses M*. Total masses were 
es timated from th e disk scale lengths and maximum rotation velocities adopting the results 
of Ivan den Bosch I (120021 1 . Our sample covers the ranges 2.0 • 10 8 M < M* < 3.7 • 10 11 M 
(median (M*) = 9.0 • 10 9 M ) in stellar mass and 2.5 • 10 9 M < M vir < 5.2 • 10 12 M (median 
(M vir ) = 1.1 -lO n M ) in total mass. 



In Fig. [51 we show the ratio between stellar and total mass as a function of redshift. If 
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Fig. 5.— The stellar-to-total mass fractions of the 110 galaxies from our sample with 
determined Vmax arid available i^-band photometry. The three large open squares show the 
median mass fraction within three redsift bins containing 36-37 galaxies each. The lack of 
a stellar mass fraction increase towards lower redshift favours a hierarchical buildup of the 
galaxies (see text for details). 
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we sub-divide our sample into three redshift bins with 36 to 37 galaxies in each bin, we find 
median stellar mass fractions of 0.068, 0.074 and 0.095 at redshifts z < 0.35, 0.35 < z < 0.6 
and z > 0.6, respectively. We hence observe a very slight decrease of the stellar mass 
fraction of spirals between z = 1 and the local universe. If late-type galaxies contained all 
their gas by redshift unity and only converted gas into stars since then, an increase of the 
stellar mass fraction would be the result. Instead of such a "monolithic" scenario, which is 
reliable to describe the evolution of massive ellipticals, the data indicate that spiral galaxies 
have accreted baryonic (and most probably also dark) mat t er in the regime < z < 1, in 
agreement with the observational findin gs of IConselice et all (120051 ) and with the expectation 



for a hierachical structure growth (e.g. iBaugh et al.ll2005l ). 



6. Conclusions 



We have used VLT/FORS spectroscopy and HST/ACS imaging to construct a sample 
of 220 field spiral galaxies up to redshift z = 1. Spatially resolved rotation curves were 
extracted and fitted with synthetic velocity fields that take into account all geometric and 
observational effects, like blurring due to the slit width and seeing influence. Using these fits, 
the maximum rotation velocity V max could be determined for 124 galaxies with an average 
look-back time of ~5 Gyr. We find that this sample is offset from the local Tully-Fisher 
Relation. The distant low-mass galaxies are more luminous at given V max than their local 
counterparts, whereas the distant high-mass spirals are compatible with the local TFR. 
Taking the magnitude limit of our sample into account, we show that the slope of local 
and distant relation would be in compliance if the TFR scatter decreased by more than a 
factor of 3 between z ~ 0.5 and z ~ 0. On the other hand, the M/L ratios indicate a 
luminosity-dependent evolution in the sense that distant low-luminosity disks have much 
lower M/L ratios than their local counterparts, while high-luminosity disks barely evolved 
in M/L over the covered redshift range. This could be interpreted as an indication of the 
"downsizing" effect, i.e. the succesive shift of star formation from high-mass to low-mass 
galaxies towards lower redshifts. In terms of the Dark Matter Halos, we find evidence for 
a hierarchical evolution, since the fraction between stellar and total mass remained roughly 
constant since z = 1. 
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